\(\int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx\) [343]

Optimal result

Integrand size = 34, antiderivative size = 517 \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=-\frac {e x}{a}+\frac {\left (1-\frac {a^2}{b^2}\right ) e x}{a}-\frac {f x^2}{2 a}+\frac {\left (1-\frac {a^2}{b^2}\right ) f x^2}{2 a}+\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}-\frac {\left (a^2-b^2\right ) (e+f x) \cos (c+d x)}{a^2 b d}-\frac {(e+f x) \cot (c+d x)}{a d}-\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {f \log (\sin (c+d x))}{a d^2}-\frac {i b f \operatorname {PolyLog}\left (2,-e^{i (c+d x)}\right )}{a^2 d^2}+\frac {i b f \operatorname {PolyLog}\left (2,e^{i (c+d x)}\right )}{a^2 d^2}-\frac {\left (a^2-b^2\right )^{3/2} f \operatorname {PolyLog}\left (2,\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d^2}+\frac {\left (a^2-b^2\right )^{3/2} f \operatorname {PolyLog}\left (2,\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d^2}+\frac {b f \sin (c+d x)}{a^2 d^2}+\frac {\left (a^2-b^2\right ) f \sin (c+d x)}{a^2 b d^2} \]

[Out]

-e*x/a+(1-a^2/b^2)*e*x/a-1/2*f*x^2/a+1/2*(1-a^2/b^2)*f*x^2/a+2*b*(f*x+e)*arctanh(exp(I*(d*x+c)))/a^2/d-b*(f*x+
e)*cos(d*x+c)/a^2/d-(a^2-b^2)*(f*x+e)*cos(d*x+c)/a^2/b/d-(f*x+e)*cot(d*x+c)/a/d+f*ln(sin(d*x+c))/a/d^2-I*(a^2-
b^2)^(3/2)*(f*x+e)*ln(1-I*b*exp(I*(d*x+c))/(a-(a^2-b^2)^(1/2)))/a^2/b^2/d-I*b*f*polylog(2,-exp(I*(d*x+c)))/a^2
/d^2+I*(a^2-b^2)^(3/2)*(f*x+e)*ln(1-I*b*exp(I*(d*x+c))/(a+(a^2-b^2)^(1/2)))/a^2/b^2/d+I*b*f*polylog(2,exp(I*(d
*x+c)))/a^2/d^2-(a^2-b^2)^(3/2)*f*polylog(2,I*b*exp(I*(d*x+c))/(a-(a^2-b^2)^(1/2)))/a^2/b^2/d^2+(a^2-b^2)^(3/2
)*f*polylog(2,I*b*exp(I*(d*x+c))/(a+(a^2-b^2)^(1/2)))/a^2/b^2/d^2+b*f*sin(d*x+c)/a^2/d^2+(a^2-b^2)*f*sin(d*x+c
)/a^2/b/d^2

Rubi [A] (verified)

Time = 1.38 (sec) , antiderivative size = 517, normalized size of antiderivative = 1.00, number of steps used = 38, number of rules used = 16, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.471, Rules used = {4639, 4493, 3391, 3801, 3556, 4490, 2713, 3377, 2717, 4268, 2317, 2438, 4621, 3404, 2296, 2221} \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {f \left (a^2-b^2\right )^{3/2} \operatorname {PolyLog}\left (2,\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d^2}+\frac {f \left (a^2-b^2\right )^{3/2} \operatorname {PolyLog}\left (2,\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d^2}+\frac {f \left (a^2-b^2\right ) \sin (c+d x)}{a^2 b d^2}-\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{\sqrt {a^2-b^2}+a}\right )}{a^2 b^2 d}-\frac {\left (a^2-b^2\right ) (e+f x) \cos (c+d x)}{a^2 b d}+\frac {e x \left (1-\frac {a^2}{b^2}\right )}{a}+\frac {f x^2 \left (1-\frac {a^2}{b^2}\right )}{2 a}-\frac {i b f \operatorname {PolyLog}\left (2,-e^{i (c+d x)}\right )}{a^2 d^2}+\frac {i b f \operatorname {PolyLog}\left (2,e^{i (c+d x)}\right )}{a^2 d^2}+\frac {b f \sin (c+d x)}{a^2 d^2}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}+\frac {f \log (\sin (c+d x))}{a d^2}-\frac {(e+f x) \cot (c+d x)}{a d}-\frac {e x}{a}-\frac {f x^2}{2 a} \]

[In]

Int[((e + f*x)*Cos[c + d*x]^2*Cot[c + d*x]^2)/(a + b*Sin[c + d*x]),x]

[Out]

-((e*x)/a) + ((1 - a^2/b^2)*e*x)/a - (f*x^2)/(2*a) + ((1 - a^2/b^2)*f*x^2)/(2*a) + (2*b*(e + f*x)*ArcTanh[E^(I
*(c + d*x))])/(a^2*d) - (b*(e + f*x)*Cos[c + d*x])/(a^2*d) - ((a^2 - b^2)*(e + f*x)*Cos[c + d*x])/(a^2*b*d) -
((e + f*x)*Cot[c + d*x])/(a*d) - (I*(a^2 - b^2)^(3/2)*(e + f*x)*Log[1 - (I*b*E^(I*(c + d*x)))/(a - Sqrt[a^2 -
b^2])])/(a^2*b^2*d) + (I*(a^2 - b^2)^(3/2)*(e + f*x)*Log[1 - (I*b*E^(I*(c + d*x)))/(a + Sqrt[a^2 - b^2])])/(a^
2*b^2*d) + (f*Log[Sin[c + d*x]])/(a*d^2) - (I*b*f*PolyLog[2, -E^(I*(c + d*x))])/(a^2*d^2) + (I*b*f*PolyLog[2,
E^(I*(c + d*x))])/(a^2*d^2) - ((a^2 - b^2)^(3/2)*f*PolyLog[2, (I*b*E^(I*(c + d*x)))/(a - Sqrt[a^2 - b^2])])/(a
^2*b^2*d^2) + ((a^2 - b^2)^(3/2)*f*PolyLog[2, (I*b*E^(I*(c + d*x)))/(a + Sqrt[a^2 - b^2])])/(a^2*b^2*d^2) + (b
*f*Sin[c + d*x])/(a^2*d^2) + ((a^2 - b^2)*f*Sin[c + d*x])/(a^2*b*d^2)

Rule 2221

Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/((a_) + (b_.)*((F_)^((g_.)*((e_.) +
(f_.)*(x_))))^(n_.)), x_Symbol] :> Simp[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x]
- Dist[d*(m/(b*f*g*n*Log[F])), Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x], x] /; FreeQ[{F,
a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]

Rule 2296

Int[((F_)^(u_)*((f_.) + (g_.)*(x_))^(m_.))/((a_.) + (b_.)*(F_)^(u_) + (c_.)*(F_)^(v_)), x_Symbol] :> With[{q =
Rt[b^2 - 4*a*c, 2]}, Dist[2*(c/q), Int[(f + g*x)^m*(F^u/(b - q + 2*c*F^u)), x], x] - Dist[2*(c/q), Int[(f + g
*x)^m*(F^u/(b + q + 2*c*F^u)), x], x]] /; FreeQ[{F, a, b, c, f, g}, x] && EqQ[v, 2*u] && LinearQ[u, x] && NeQ[
b^2 - 4*a*c, 0] && IGtQ[m, 0]

Rule 2317

Int[Log[(a_) + (b_.)*((F_)^((e_.)*((c_.) + (d_.)*(x_))))^(n_.)], x_Symbol] :> Dist[1/(d*e*n*Log[F]), Subst[Int
[Log[a + b*x]/x, x], x, (F^(e*(c + d*x)))^n], x] /; FreeQ[{F, a, b, c, d, e, n}, x] && GtQ[a, 0]

Rule 2438

Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2, (-c)*e*x^n]/n, x] /; FreeQ[{c, d,
e, n}, x] && EqQ[c*d, 1]

Rule 2713

Int[sin[(c_.) + (d_.)*(x_)]^(n_), x_Symbol] :> Dist[-d^(-1), Subst[Int[Expand[(1 - x^2)^((n - 1)/2), x], x], x
, Cos[c + d*x]], x] /; FreeQ[{c, d}, x] && IGtQ[(n - 1)/2, 0]

Rule 2717

Int[sin[Pi/2 + (c_.) + (d_.)*(x_)], x_Symbol] :> Simp[Sin[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3377

Int[((c_.) + (d_.)*(x_))^(m_.)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[(-(c + d*x)^m)*(Cos[e + f*x]/f), x]
+ Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && GtQ[m, 0]

Rule 3391

Int[((c_.) + (d_.)*(x_))*((b_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[d*((b*Sin[e + f*x])^n/(f^2*n^
2)), x] + (Dist[b^2*((n - 1)/n), Int[(c + d*x)*(b*Sin[e + f*x])^(n - 2), x], x] - Simp[b*(c + d*x)*Cos[e + f*x
]*((b*Sin[e + f*x])^(n - 1)/(f*n)), x]) /; FreeQ[{b, c, d, e, f}, x] && GtQ[n, 1]

Rule 3404

Int[((c_.) + (d_.)*(x_))^(m_.)/((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Dist[2, Int[(c + d*x)^m*(E
^(I*(e + f*x))/(I*b + 2*a*E^(I*(e + f*x)) - I*b*E^(2*I*(e + f*x)))), x], x] /; FreeQ[{a, b, c, d, e, f}, x] &&
NeQ[a^2 - b^2, 0] && IGtQ[m, 0]

Rule 3556

Int[tan[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Log[RemoveContent[Cos[c + d*x], x]]/d, x] /; FreeQ[{c, d}, x]

Rule 3801

Int[((c_.) + (d_.)*(x_))^(m_.)*((b_.)*tan[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[b*(c + d*x)^m*((b*Tan[e
+ f*x])^(n - 1)/(f*(n - 1))), x] + (-Dist[b*d*(m/(f*(n - 1))), Int[(c + d*x)^(m - 1)*(b*Tan[e + f*x])^(n - 1)
, x], x] - Dist[b^2, Int[(c + d*x)^m*(b*Tan[e + f*x])^(n - 2), x], x]) /; FreeQ[{b, c, d, e, f}, x] && GtQ[n,
1] && GtQ[m, 0]

Rule 4268

Int[csc[(e_.) + (f_.)*(x_)]*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> Simp[-2*(c + d*x)^m*(ArcTanh[E^(I*(e + f*
x))]/f), x] + (-Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Log[1 - E^(I*(e + f*x))], x], x] + Dist[d*(m/f), Int[(c +
d*x)^(m - 1)*Log[1 + E^(I*(e + f*x))], x], x]) /; FreeQ[{c, d, e, f}, x] && IGtQ[m, 0]

Rule 4490

Int[Cos[(a_.) + (b_.)*(x_)]^(n_.)*((c_.) + (d_.)*(x_))^(m_.)*Sin[(a_.) + (b_.)*(x_)], x_Symbol] :> Simp[(-(c +
d*x)^m)*(Cos[a + b*x]^(n + 1)/(b*(n + 1))), x] + Dist[d*(m/(b*(n + 1))), Int[(c + d*x)^(m - 1)*Cos[a + b*x]^(
n + 1), x], x] /; FreeQ[{a, b, c, d, n}, x] && IGtQ[m, 0] && NeQ[n, -1]

Rule 4493

Int[Cos[(a_.) + (b_.)*(x_)]^(n_.)*Cot[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> -Int[
(c + d*x)^m*Cos[a + b*x]^n*Cot[a + b*x]^(p - 2), x] + Int[(c + d*x)^m*Cos[a + b*x]^(n - 2)*Cot[a + b*x]^p, x]
/; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0] && IGtQ[p, 0]

Rule 4621

Int[(Cos[(c_.) + (d_.)*(x_)]^(n_)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sin[(c_.) + (d_.)*(x_)]), x_Symbol
] :> Dist[a/b^2, Int[(e + f*x)^m*Cos[c + d*x]^(n - 2), x], x] + (-Dist[1/b, Int[(e + f*x)^m*Cos[c + d*x]^(n -
2)*Sin[c + d*x], x], x] - Dist[(a^2 - b^2)/b^2, Int[(e + f*x)^m*(Cos[c + d*x]^(n - 2)/(a + b*Sin[c + d*x])), x
], x]) /; FreeQ[{a, b, c, d, e, f}, x] && IGtQ[n, 1] && NeQ[a^2 - b^2, 0] && IGtQ[m, 0]

Rule 4639

Int[(Cos[(c_.) + (d_.)*(x_)]^(p_.)*Cot[(c_.) + (d_.)*(x_)]^(n_.)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sin
[(c_.) + (d_.)*(x_)]), x_Symbol] :> Dist[1/a, Int[(e + f*x)^m*Cos[c + d*x]^p*Cot[c + d*x]^n, x], x] - Dist[b/a
, Int[(e + f*x)^m*Cos[c + d*x]^(p + 1)*(Cot[c + d*x]^(n - 1)/(a + b*Sin[c + d*x])), x], x] /; FreeQ[{a, b, c,
d, e, f}, x] && IGtQ[m, 0] && IGtQ[n, 0] && IGtQ[p, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\int (e+f x) \cos ^2(c+d x) \cot ^2(c+d x) \, dx}{a}-\frac {b \int \frac {(e+f x) \cos ^3(c+d x) \cot (c+d x)}{a+b \sin (c+d x)} \, dx}{a} \\ & = -\frac {\int (e+f x) \cos ^2(c+d x) \, dx}{a}+\frac {\int (e+f x) \cot ^2(c+d x) \, dx}{a}-\frac {b \int (e+f x) \cos ^3(c+d x) \cot (c+d x) \, dx}{a^2}+\frac {b^2 \int \frac {(e+f x) \cos ^4(c+d x)}{a+b \sin (c+d x)} \, dx}{a^2} \\ & = -\frac {f \cos ^2(c+d x)}{4 a d^2}-\frac {(e+f x) \cot (c+d x)}{a d}-\frac {(e+f x) \cos (c+d x) \sin (c+d x)}{2 a d}-\frac {\int (e+f x) \, dx}{2 a}-\frac {\int (e+f x) \, dx}{a}+\frac {\int (e+f x) \cos ^2(c+d x) \, dx}{a}-\frac {b \int (e+f x) \cos (c+d x) \cot (c+d x) \, dx}{a^2}-\left (1-\frac {b^2}{a^2}\right ) \int \frac {(e+f x) \cos ^2(c+d x)}{a+b \sin (c+d x)} \, dx+\frac {f \int \cot (c+d x) \, dx}{a d} \\ & = -\frac {3 e x}{2 a}-\frac {3 f x^2}{4 a}-\frac {(e+f x) \cot (c+d x)}{a d}+\frac {f \log (\sin (c+d x))}{a d^2}+\frac {\int (e+f x) \, dx}{2 a}-\frac {b \int (e+f x) \csc (c+d x) \, dx}{a^2}+\frac {b \int (e+f x) \sin (c+d x) \, dx}{a^2}-\frac {\left (a \left (1-\frac {b^2}{a^2}\right )\right ) \int (e+f x) \, dx}{b^2}-\frac {\left (-1+\frac {b^2}{a^2}\right ) \int (e+f x) \sin (c+d x) \, dx}{b}-\frac {\left (\left (a^2-b^2\right ) \left (-1+\frac {b^2}{a^2}\right )\right ) \int \frac {e+f x}{a+b \sin (c+d x)} \, dx}{b^2} \\ & = -\frac {e x}{a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) e x}{b^2}-\frac {f x^2}{2 a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) f x^2}{2 b^2}+\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}-\frac {\left (1-\frac {b^2}{a^2}\right ) (e+f x) \cos (c+d x)}{b d}-\frac {(e+f x) \cot (c+d x)}{a d}+\frac {f \log (\sin (c+d x))}{a d^2}+\frac {\left (2 \left (a^2-b^2\right )^2\right ) \int \frac {e^{i (c+d x)} (e+f x)}{i b+2 a e^{i (c+d x)}-i b e^{2 i (c+d x)}} \, dx}{a^2 b^2}+\frac {(b f) \int \cos (c+d x) \, dx}{a^2 d}+\frac {(b f) \int \log \left (1-e^{i (c+d x)}\right ) \, dx}{a^2 d}-\frac {(b f) \int \log \left (1+e^{i (c+d x)}\right ) \, dx}{a^2 d}+\frac {\left (\left (1-\frac {b^2}{a^2}\right ) f\right ) \int \cos (c+d x) \, dx}{b d} \\ & = -\frac {e x}{a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) e x}{b^2}-\frac {f x^2}{2 a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) f x^2}{2 b^2}+\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}-\frac {\left (1-\frac {b^2}{a^2}\right ) (e+f x) \cos (c+d x)}{b d}-\frac {(e+f x) \cot (c+d x)}{a d}+\frac {f \log (\sin (c+d x))}{a d^2}+\frac {b f \sin (c+d x)}{a^2 d^2}+\frac {\left (1-\frac {b^2}{a^2}\right ) f \sin (c+d x)}{b d^2}-\frac {\left (2 i \left (a^2-b^2\right )^{3/2}\right ) \int \frac {e^{i (c+d x)} (e+f x)}{2 a-2 \sqrt {a^2-b^2}-2 i b e^{i (c+d x)}} \, dx}{a^2 b}+\frac {\left (2 i \left (a^2-b^2\right )^{3/2}\right ) \int \frac {e^{i (c+d x)} (e+f x)}{2 a+2 \sqrt {a^2-b^2}-2 i b e^{i (c+d x)}} \, dx}{a^2 b}-\frac {(i b f) \text {Subst}\left (\int \frac {\log (1-x)}{x} \, dx,x,e^{i (c+d x)}\right )}{a^2 d^2}+\frac {(i b f) \text {Subst}\left (\int \frac {\log (1+x)}{x} \, dx,x,e^{i (c+d x)}\right )}{a^2 d^2} \\ & = -\frac {e x}{a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) e x}{b^2}-\frac {f x^2}{2 a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) f x^2}{2 b^2}+\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}-\frac {\left (1-\frac {b^2}{a^2}\right ) (e+f x) \cos (c+d x)}{b d}-\frac {(e+f x) \cot (c+d x)}{a d}-\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {f \log (\sin (c+d x))}{a d^2}-\frac {i b f \operatorname {PolyLog}\left (2,-e^{i (c+d x)}\right )}{a^2 d^2}+\frac {i b f \operatorname {PolyLog}\left (2,e^{i (c+d x)}\right )}{a^2 d^2}+\frac {b f \sin (c+d x)}{a^2 d^2}+\frac {\left (1-\frac {b^2}{a^2}\right ) f \sin (c+d x)}{b d^2}+\frac {\left (i \left (a^2-b^2\right )^{3/2} f\right ) \int \log \left (1-\frac {2 i b e^{i (c+d x)}}{2 a-2 \sqrt {a^2-b^2}}\right ) \, dx}{a^2 b^2 d}-\frac {\left (i \left (a^2-b^2\right )^{3/2} f\right ) \int \log \left (1-\frac {2 i b e^{i (c+d x)}}{2 a+2 \sqrt {a^2-b^2}}\right ) \, dx}{a^2 b^2 d} \\ & = -\frac {e x}{a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) e x}{b^2}-\frac {f x^2}{2 a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) f x^2}{2 b^2}+\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}-\frac {\left (1-\frac {b^2}{a^2}\right ) (e+f x) \cos (c+d x)}{b d}-\frac {(e+f x) \cot (c+d x)}{a d}-\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {f \log (\sin (c+d x))}{a d^2}-\frac {i b f \operatorname {PolyLog}\left (2,-e^{i (c+d x)}\right )}{a^2 d^2}+\frac {i b f \operatorname {PolyLog}\left (2,e^{i (c+d x)}\right )}{a^2 d^2}+\frac {b f \sin (c+d x)}{a^2 d^2}+\frac {\left (1-\frac {b^2}{a^2}\right ) f \sin (c+d x)}{b d^2}+\frac {\left (\left (a^2-b^2\right )^{3/2} f\right ) \text {Subst}\left (\int \frac {\log \left (1-\frac {2 i b x}{2 a-2 \sqrt {a^2-b^2}}\right )}{x} \, dx,x,e^{i (c+d x)}\right )}{a^2 b^2 d^2}-\frac {\left (\left (a^2-b^2\right )^{3/2} f\right ) \text {Subst}\left (\int \frac {\log \left (1-\frac {2 i b x}{2 a+2 \sqrt {a^2-b^2}}\right )}{x} \, dx,x,e^{i (c+d x)}\right )}{a^2 b^2 d^2} \\ & = -\frac {e x}{a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) e x}{b^2}-\frac {f x^2}{2 a}-\frac {a \left (1-\frac {b^2}{a^2}\right ) f x^2}{2 b^2}+\frac {2 b (e+f x) \text {arctanh}\left (e^{i (c+d x)}\right )}{a^2 d}-\frac {b (e+f x) \cos (c+d x)}{a^2 d}-\frac {\left (1-\frac {b^2}{a^2}\right ) (e+f x) \cos (c+d x)}{b d}-\frac {(e+f x) \cot (c+d x)}{a d}-\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {i \left (a^2-b^2\right )^{3/2} (e+f x) \log \left (1-\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d}+\frac {f \log (\sin (c+d x))}{a d^2}-\frac {i b f \operatorname {PolyLog}\left (2,-e^{i (c+d x)}\right )}{a^2 d^2}+\frac {i b f \operatorname {PolyLog}\left (2,e^{i (c+d x)}\right )}{a^2 d^2}-\frac {\left (a^2-b^2\right )^{3/2} f \operatorname {PolyLog}\left (2,\frac {i b e^{i (c+d x)}}{a-\sqrt {a^2-b^2}}\right )}{a^2 b^2 d^2}+\frac {\left (a^2-b^2\right )^{3/2} f \operatorname {PolyLog}\left (2,\frac {i b e^{i (c+d x)}}{a+\sqrt {a^2-b^2}}\right )}{a^2 b^2 d^2}+\frac {b f \sin (c+d x)}{a^2 d^2}+\frac {\left (1-\frac {b^2}{a^2}\right ) f \sin (c+d x)}{b d^2} \\ \end{align*}

Mathematica [B] (warning: unable to verify)

Both result and optimal contain complex but leaf count is larger than twice the leaf count of optimal. \(1091\) vs. \(2(517)=1034\).

Time = 12.25 (sec) , antiderivative size = 1091, normalized size of antiderivative = 2.11 \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=-\frac {a (c+d x) (2 d e-2 c f+f (c+d x))}{2 b^2 d^2}-\frac {(d e-c f+f (c+d x)) \cos (c+d x)}{b d^2}+\frac {\left (-d e \cos \left (\frac {1}{2} (c+d x)\right )+c f \cos \left (\frac {1}{2} (c+d x)\right )-f (c+d x) \cos \left (\frac {1}{2} (c+d x)\right )\right ) \csc \left (\frac {1}{2} (c+d x)\right )}{2 a d^2}-\frac {b e \log \left (\tan \left (\frac {1}{2} (c+d x)\right )\right )}{a^2 d}+\frac {b c f \log \left (\tan \left (\frac {1}{2} (c+d x)\right )\right )}{a^2 d^2}+\frac {f (\log (\cos (c+d x))+\log (\tan (c+d x)))}{a d^2}-\frac {b f \left ((c+d x) \left (\log \left (1-e^{i (c+d x)}\right )-\log \left (1+e^{i (c+d x)}\right )\right )+i \left (\operatorname {PolyLog}\left (2,-e^{i (c+d x)}\right )-\operatorname {PolyLog}\left (2,e^{i (c+d x)}\right )\right )\right )}{a^2 d^2}+\frac {\left (a^2-b^2\right )^2 (d e+d f x) \left (\frac {2 (d e-c f) \arctan \left (\frac {b+a \tan \left (\frac {1}{2} (c+d x)\right )}{\sqrt {a^2-b^2}}\right )}{\sqrt {a^2-b^2}}-\frac {i f \log \left (1+i \tan \left (\frac {1}{2} (c+d x)\right )\right ) \log \left (\frac {b-\sqrt {-a^2+b^2}+a \tan \left (\frac {1}{2} (c+d x)\right )}{i a+b-\sqrt {-a^2+b^2}}\right )}{\sqrt {-a^2+b^2}}+\frac {i f \log \left (1-i \tan \left (\frac {1}{2} (c+d x)\right )\right ) \log \left (-\frac {b-\sqrt {-a^2+b^2}+a \tan \left (\frac {1}{2} (c+d x)\right )}{i a-b+\sqrt {-a^2+b^2}}\right )}{\sqrt {-a^2+b^2}}-\frac {i f \log \left (1-i \tan \left (\frac {1}{2} (c+d x)\right )\right ) \log \left (\frac {b+\sqrt {-a^2+b^2}+a \tan \left (\frac {1}{2} (c+d x)\right )}{-i a+b+\sqrt {-a^2+b^2}}\right )}{\sqrt {-a^2+b^2}}+\frac {i f \log \left (1+i \tan \left (\frac {1}{2} (c+d x)\right )\right ) \log \left (\frac {b+\sqrt {-a^2+b^2}+a \tan \left (\frac {1}{2} (c+d x)\right )}{i a+b+\sqrt {-a^2+b^2}}\right )}{\sqrt {-a^2+b^2}}-\frac {i f \operatorname {PolyLog}\left (2,\frac {a \left (1-i \tan \left (\frac {1}{2} (c+d x)\right )\right )}{a+i \left (b+\sqrt {-a^2+b^2}\right )}\right )}{\sqrt {-a^2+b^2}}+\frac {i f \operatorname {PolyLog}\left (2,\frac {a \left (1+i \tan \left (\frac {1}{2} (c+d x)\right )\right )}{a-i \left (b+\sqrt {-a^2+b^2}\right )}\right )}{\sqrt {-a^2+b^2}}+\frac {i f \operatorname {PolyLog}\left (2,\frac {a \left (i+\tan \left (\frac {1}{2} (c+d x)\right )\right )}{i a-b+\sqrt {-a^2+b^2}}\right )}{\sqrt {-a^2+b^2}}-\frac {i f \operatorname {PolyLog}\left (2,\frac {a+i a \tan \left (\frac {1}{2} (c+d x)\right )}{a+i \left (-b+\sqrt {-a^2+b^2}\right )}\right )}{\sqrt {-a^2+b^2}}\right )}{a^2 b^2 d^2 \left (d e-c f+i f \log \left (1-i \tan \left (\frac {1}{2} (c+d x)\right )\right )-i f \log \left (1+i \tan \left (\frac {1}{2} (c+d x)\right )\right )\right )}+\frac {\sec \left (\frac {1}{2} (c+d x)\right ) \left (d e \sin \left (\frac {1}{2} (c+d x)\right )-c f \sin \left (\frac {1}{2} (c+d x)\right )+f (c+d x) \sin \left (\frac {1}{2} (c+d x)\right )\right )}{2 a d^2}+\frac {f \sin (c+d x)}{b d^2} \]

[In]

Integrate[((e + f*x)*Cos[c + d*x]^2*Cot[c + d*x]^2)/(a + b*Sin[c + d*x]),x]

[Out]

-1/2*(a*(c + d*x)*(2*d*e - 2*c*f + f*(c + d*x)))/(b^2*d^2) - ((d*e - c*f + f*(c + d*x))*Cos[c + d*x])/(b*d^2)
+ ((-(d*e*Cos[(c + d*x)/2]) + c*f*Cos[(c + d*x)/2] - f*(c + d*x)*Cos[(c + d*x)/2])*Csc[(c + d*x)/2])/(2*a*d^2)
- (b*e*Log[Tan[(c + d*x)/2]])/(a^2*d) + (b*c*f*Log[Tan[(c + d*x)/2]])/(a^2*d^2) + (f*(Log[Cos[c + d*x]] + Log
[Tan[c + d*x]]))/(a*d^2) - (b*f*((c + d*x)*(Log[1 - E^(I*(c + d*x))] - Log[1 + E^(I*(c + d*x))]) + I*(PolyLog[
2, -E^(I*(c + d*x))] - PolyLog[2, E^(I*(c + d*x))])))/(a^2*d^2) + ((a^2 - b^2)^2*(d*e + d*f*x)*((2*(d*e - c*f)
*ArcTan[(b + a*Tan[(c + d*x)/2])/Sqrt[a^2 - b^2]])/Sqrt[a^2 - b^2] - (I*f*Log[1 + I*Tan[(c + d*x)/2]]*Log[(b -
Sqrt[-a^2 + b^2] + a*Tan[(c + d*x)/2])/(I*a + b - Sqrt[-a^2 + b^2])])/Sqrt[-a^2 + b^2] + (I*f*Log[1 - I*Tan[(
c + d*x)/2]]*Log[-((b - Sqrt[-a^2 + b^2] + a*Tan[(c + d*x)/2])/(I*a - b + Sqrt[-a^2 + b^2]))])/Sqrt[-a^2 + b^2
] - (I*f*Log[1 - I*Tan[(c + d*x)/2]]*Log[(b + Sqrt[-a^2 + b^2] + a*Tan[(c + d*x)/2])/((-I)*a + b + Sqrt[-a^2 +
b^2])])/Sqrt[-a^2 + b^2] + (I*f*Log[1 + I*Tan[(c + d*x)/2]]*Log[(b + Sqrt[-a^2 + b^2] + a*Tan[(c + d*x)/2])/(
I*a + b + Sqrt[-a^2 + b^2])])/Sqrt[-a^2 + b^2] - (I*f*PolyLog[2, (a*(1 - I*Tan[(c + d*x)/2]))/(a + I*(b + Sqrt
[-a^2 + b^2]))])/Sqrt[-a^2 + b^2] + (I*f*PolyLog[2, (a*(1 + I*Tan[(c + d*x)/2]))/(a - I*(b + Sqrt[-a^2 + b^2])
)])/Sqrt[-a^2 + b^2] + (I*f*PolyLog[2, (a*(I + Tan[(c + d*x)/2]))/(I*a - b + Sqrt[-a^2 + b^2])])/Sqrt[-a^2 + b
^2] - (I*f*PolyLog[2, (a + I*a*Tan[(c + d*x)/2])/(a + I*(-b + Sqrt[-a^2 + b^2]))])/Sqrt[-a^2 + b^2]))/(a^2*b^2
*d^2*(d*e - c*f + I*f*Log[1 - I*Tan[(c + d*x)/2]] - I*f*Log[1 + I*Tan[(c + d*x)/2]])) + (Sec[(c + d*x)/2]*(d*e
*Sin[(c + d*x)/2] - c*f*Sin[(c + d*x)/2] + f*(c + d*x)*Sin[(c + d*x)/2]))/(2*a*d^2) + (f*Sin[c + d*x])/(b*d^2)

Maple [B] (verified)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 5309 vs. \(2 (475 ) = 950\).

Time = 2.82 (sec) , antiderivative size = 5310, normalized size of antiderivative = 10.27

method	result	size
risch	\(\text {Expression too large to display}\)	\(5310\)

[In]

int((f*x+e)*cos(d*x+c)^2*cot(d*x+c)^2/(a+b*sin(d*x+c)),x,method=_RETURNVERBOSE)

[Out]

result too large to display

Fricas [B] (verification not implemented)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 1751 vs. \(2 (462) = 924\).

Time = 0.55 (sec) , antiderivative size = 1751, normalized size of antiderivative = 3.39 \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=\text {Too large to display} \]

[In]

integrate((f*x+e)*cos(d*x+c)^2*cot(d*x+c)^2/(a+b*sin(d*x+c)),x, algorithm="fricas")

[Out]

-1/2*(2*a^2*b*f*cos(d*x + c)^2 - I*b^3*f*dilog(cos(d*x + c) + I*sin(d*x + c))*sin(d*x + c) + I*b^3*f*dilog(cos
(d*x + c) - I*sin(d*x + c))*sin(d*x + c) - I*b^3*f*dilog(-cos(d*x + c) + I*sin(d*x + c))*sin(d*x + c) + I*b^3*
f*dilog(-cos(d*x + c) - I*sin(d*x + c))*sin(d*x + c) - I*(a^2*b - b^3)*f*sqrt(-(a^2 - b^2)/b^2)*dilog((I*a*cos
(d*x + c) - a*sin(d*x + c) + (b*cos(d*x + c) + I*b*sin(d*x + c))*sqrt(-(a^2 - b^2)/b^2) - b)/b + 1)*sin(d*x +
c) + I*(a^2*b - b^3)*f*sqrt(-(a^2 - b^2)/b^2)*dilog((I*a*cos(d*x + c) - a*sin(d*x + c) - (b*cos(d*x + c) + I*b
*sin(d*x + c))*sqrt(-(a^2 - b^2)/b^2) - b)/b + 1)*sin(d*x + c) + I*(a^2*b - b^3)*f*sqrt(-(a^2 - b^2)/b^2)*dilo
g((-I*a*cos(d*x + c) - a*sin(d*x + c) + (b*cos(d*x + c) - I*b*sin(d*x + c))*sqrt(-(a^2 - b^2)/b^2) - b)/b + 1)
*sin(d*x + c) - I*(a^2*b - b^3)*f*sqrt(-(a^2 - b^2)/b^2)*dilog((-I*a*cos(d*x + c) - a*sin(d*x + c) - (b*cos(d*
x + c) - I*b*sin(d*x + c))*sqrt(-(a^2 - b^2)/b^2) - b)/b + 1)*sin(d*x + c) - 2*a^2*b*f - ((a^2*b - b^3)*d*e -
(a^2*b - b^3)*c*f)*sqrt(-(a^2 - b^2)/b^2)*log(2*b*cos(d*x + c) + 2*I*b*sin(d*x + c) + 2*b*sqrt(-(a^2 - b^2)/b^
2) + 2*I*a)*sin(d*x + c) - ((a^2*b - b^3)*d*e - (a^2*b - b^3)*c*f)*sqrt(-(a^2 - b^2)/b^2)*log(2*b*cos(d*x + c)
- 2*I*b*sin(d*x + c) + 2*b*sqrt(-(a^2 - b^2)/b^2) - 2*I*a)*sin(d*x + c) + ((a^2*b - b^3)*d*e - (a^2*b - b^3)*
c*f)*sqrt(-(a^2 - b^2)/b^2)*log(-2*b*cos(d*x + c) + 2*I*b*sin(d*x + c) + 2*b*sqrt(-(a^2 - b^2)/b^2) + 2*I*a)*s
in(d*x + c) + ((a^2*b - b^3)*d*e - (a^2*b - b^3)*c*f)*sqrt(-(a^2 - b^2)/b^2)*log(-2*b*cos(d*x + c) - 2*I*b*sin
(d*x + c) + 2*b*sqrt(-(a^2 - b^2)/b^2) - 2*I*a)*sin(d*x + c) + ((a^2*b - b^3)*d*f*x + (a^2*b - b^3)*c*f)*sqrt(
-(a^2 - b^2)/b^2)*log(-(I*a*cos(d*x + c) - a*sin(d*x + c) + (b*cos(d*x + c) + I*b*sin(d*x + c))*sqrt(-(a^2 - b
^2)/b^2) - b)/b)*sin(d*x + c) - ((a^2*b - b^3)*d*f*x + (a^2*b - b^3)*c*f)*sqrt(-(a^2 - b^2)/b^2)*log(-(I*a*cos
(d*x + c) - a*sin(d*x + c) - (b*cos(d*x + c) + I*b*sin(d*x + c))*sqrt(-(a^2 - b^2)/b^2) - b)/b)*sin(d*x + c) +
((a^2*b - b^3)*d*f*x + (a^2*b - b^3)*c*f)*sqrt(-(a^2 - b^2)/b^2)*log(-(-I*a*cos(d*x + c) - a*sin(d*x + c) + (
b*cos(d*x + c) - I*b*sin(d*x + c))*sqrt(-(a^2 - b^2)/b^2) - b)/b)*sin(d*x + c) - ((a^2*b - b^3)*d*f*x + (a^2*b
- b^3)*c*f)*sqrt(-(a^2 - b^2)/b^2)*log(-(-I*a*cos(d*x + c) - a*sin(d*x + c) - (b*cos(d*x + c) - I*b*sin(d*x +
c))*sqrt(-(a^2 - b^2)/b^2) - b)/b)*sin(d*x + c) - (b^3*d*f*x + b^3*d*e + a*b^2*f)*log(cos(d*x + c) + I*sin(d*
x + c) + 1)*sin(d*x + c) - (b^3*d*f*x + b^3*d*e + a*b^2*f)*log(cos(d*x + c) - I*sin(d*x + c) + 1)*sin(d*x + c)
+ (b^3*d*e - (b^3*c + a*b^2)*f)*log(-1/2*cos(d*x + c) + 1/2*I*sin(d*x + c) + 1/2)*sin(d*x + c) + (b^3*d*e - (
b^3*c + a*b^2)*f)*log(-1/2*cos(d*x + c) - 1/2*I*sin(d*x + c) + 1/2)*sin(d*x + c) + (b^3*d*f*x + b^3*c*f)*log(-
cos(d*x + c) + I*sin(d*x + c) + 1)*sin(d*x + c) + (b^3*d*f*x + b^3*c*f)*log(-cos(d*x + c) - I*sin(d*x + c) + 1
)*sin(d*x + c) + 2*(a*b^2*d*f*x + a*b^2*d*e)*cos(d*x + c) + (a^3*d^2*f*x^2 + 2*a^3*d^2*e*x + 2*(a^2*b*d*f*x +
a^2*b*d*e)*cos(d*x + c))*sin(d*x + c))/(a^2*b^2*d^2*sin(d*x + c))

Sympy [F]

\[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=\int \frac {\left (e + f x\right ) \cos ^{2}{\left (c + d x \right )} \cot ^{2}{\left (c + d x \right )}}{a + b \sin {\left (c + d x \right )}}\, dx \]

[In]

integrate((f*x+e)*cos(d*x+c)**2*cot(d*x+c)**2/(a+b*sin(d*x+c)),x)

[Out]

Integral((e + f*x)*cos(c + d*x)**2*cot(c + d*x)**2/(a + b*sin(c + d*x)), x)

Maxima [F(-2)]

Exception generated. \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=\text {Exception raised: ValueError} \]

[In]

integrate((f*x+e)*cos(d*x+c)^2*cot(d*x+c)^2/(a+b*sin(d*x+c)),x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(4*b^2-4*a^2>0)', see `assume?`
for more de

Giac [F(-1)]

Timed out. \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=\text {Timed out} \]

[In]

integrate((f*x+e)*cos(d*x+c)^2*cot(d*x+c)^2/(a+b*sin(d*x+c)),x, algorithm="giac")

[Out]

Timed out

Mupad [F(-1)]

Timed out. \[ \int \frac {(e+f x) \cos ^2(c+d x) \cot ^2(c+d x)}{a+b \sin (c+d x)} \, dx=\text {Hanged} \]

[In]

int((cos(c + d*x)^2*cot(c + d*x)^2*(e + f*x))/(a + b*sin(c + d*x)),x)

[Out]

\text{Hanged}